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Abstract 

In this paper, the transition of synchronizing path of delay-coupled chaotic 
oscillators in a scale-free network is highlighted. Mainly, through the critical 
transmission delay makes chaotic oscillators be coupled on the edge of sta¬ 
bility, we find that the transition of synchronizing path is abnormal , which is 
characterized by the following evidences: (a) synchronization process starts 
with low-degree rather than high-degree ones; (b) the high-degree nodes don’t 
undertake the role of hub; (c) the synchronized subnetworks show a poor 
small-world property as a result of hubs absence; (d) the clustering syn¬ 
chronization behavior emerges even community structure is absent in the 
scale-free network. This abnormal synchronizing path suggests that the di¬ 
verse synchronization behaviors occur in the same topology, which implies 
that the relationship between dynamics and structure of network is much 
more complicated than the common sense that the structure is the founda¬ 
tion of dynamics. Moreover, it also reveals the potential connection from the 
transition of synchronization behavior to disorder in real complex networks, 
e.g. Alzheimer disease 
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1. Introduction 


Synchronization, which is referred to as the uniform activities of inter¬ 
acting units, is ubiquitous in nature and plays an important role in physics, 
biology, sociology, technology and neural systems 0,®. Generally, real-world 
systems are often represented by networks ,of which the nodes are units and 
edges indicate their interactions 0, B S jb]. A networked system consist¬ 
ing of N coupled oscillators (i.e., nodes) thus can be simply described by 
x I = F( Xi )-cE,',,H( 

Xj), where x* describes the the m-dimensional state of 
node i (i = 1, 2,..., N), F indicates nodes’ dynamics, c and F are the coupling 
strength and function, and L = (l t j)wxN is the Laplacian matrix of network. 
It reaches synchronization if all oscillators satisfy x^f) = x 2 (f) = ... = Xjv(t). 

The synchronizability of network usually emphasizes whether its synchro¬ 
nization is achievable. Pecora and Carroll 0] varied coupling strength and 
firstly unraveled that the synchronizability of network could be evaluated by 
the eigenratio X ma x/^ 2 , where \ max and A 2 were the maximum and mini¬ 
mum non-trivial eigenvalue of L. The smaller the eigenratio is, the stronger 
synchronizability the network has, and vise versa. Following this pioneer¬ 
ing work, many researches have comprehensive investigated how to improve 
synchronizability of network via edge betweenness 01, to pology modification 
[a, 3, optimization 11, 12, 13||, adaptive evolution pB, and so on. Besides 
the network structure, Li and Chen 0] found that transmission delay be¬ 
tween oscillators also influenced synchronizability of network deeply due to 
the signals traveling in limited speed, such as neural oscillators |16| . 

Moreover, the temporal characteristics shown in synchronization process 
can provide more information of the structure and function of system. For 
examples, the temporal correlations between interacting neural oscillators 


represent the clustering structure of coupling topology [17|, [18 


and clus¬ 
tering synchronization has been applied in community detection 19), [20] and 
data clustering (21,22]. It is also worthy to discuss one of the most attractive 
temporal characteristics, synchronizing path that describes how the synchro¬ 
nized state diffuses to whole network. Gomez-Gardenes et al 0 studied 
the synchronization path based on Kuramoto model, and showed it behave 
difference from Barabasi-Albert scale-free (BA) network to Erdos-Renyi ran¬ 
dom (ER) network. Especially, the high-degree nodes in BA network were 
more likely to be synchronized into the central cluster than low-degree ones. 
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And, Stout et al 24j proved that the synchronizing path observed in Ref. 

23] was robust to increasing oscillator number and changing the distribution 
function of the Kuramoto oscillators’ natural frequencies. Recently, Zhou 
et.al 25| studied the influence of the mixing pattern on synchronizing path 
and found that the disassortative BA networks synchronize via a merging 
process rather than the typical growth process on BA network. 

In this paper, we mainly investigate how the transmission delay affects 
the synchronization of chaotic oscillators. Based on a BA network consisting 
of Rossler oscillators, two distinct types of synchronizing paths are found via 
controlling transmission delay. Especially, when transmission delay achieves 
a critical value, the synchronizing path shows abnormal evidences that (a) 
these low-degree nodes rather than high-degree ones firstly synchronized, (b) 
the high-degree nodes no longer act as hubs in synchronization process be¬ 
cause they join the largest synchronized cluster very late, (c) the largest syn¬ 
chronized cluster show poor small-world property as a result of absent hubs, 
(d) the clustering synchronization behaviors appears even the BA network is 
in absence of community. These features make the abnormal synchronizing 
path on the edge of stability (i.e., critical transmission delay) obviously dif¬ 
fer from that on the stable region, and also reveal the disagreement between 
structure and dynamics of network. Meanwhile, this transition of synchroniz¬ 
ing path suggests that the same topology can display diverse synchronization 
behaviors, which may provide a deep insight to understand the Alzheimer 
disease. 


2. Synchronization of Delay-coupled Chaotic Oscillators 

The network consisting of N delay-coupled oscillators is represented by 

*»(*) = F M*)] - - r)], (1) 

where r is the transmission delay, c is the coupling strength, and Uj is the 
entry of L at the i-th line and j'-th column. The coupling function H(x) is 
set in the linear form as H(x) = x. As proved in Ref. |26|, there exists a 
critical transmission delay r c that the synchronization is stable (unstable) if 
t < r c (r > r c ). The synchronization error W{t) is adopted to measure the 
degree of synchronization, 

W(t) = - (x,)| 2 , (2) 
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where (...) denotes averaging over all nodes. Then the logarithmic synchro¬ 
nizing speed is defined as 


m - 0 ) 

Positive synchronizing speed means that the W(t) vanishes and the synchro¬ 
nization is stable. Thus, the critical transmission delay can be found by 
checking the sign of synchronizing speeds. 

The Rbssler oscillator is chosen for the numeric estimation of synchronous 
behaviors, which reads as 


{ i'il %i2 

Xi 2 = Xu + 0.2x i2 , (4) 

x i3 = 0.2 + I j3 (xji - 12). 


We use the Barabasi-Albert (BA) model with node number N = 1000 and 
mean degree k = 10 as the coupling topology and set coupling strength as 
c = 3. The synchronization processes are carried out by fourth-order Ronge- 
Kutta method with fixed step 5t = 1.0 x 10~ 4 . 

The synchronizing speeds with varying transmission delays are obtained 
by the first-order polynomial fitting of ln[W(t)\ after transient time of 10 4 

It can be seen that the critical transmis- 
Then the synchronization with r c = 4.2 x ICO 3 


steps, which is shown in Fig. 1(a) 
sion delay is around 4.2 x ICR 3 . 


act as an example of collective dynamics on the edge of stability and the one 
with r = 2.0 x 10~ 3 as that in the stable region to be compared with. Both 


of the temporal evolutions of W(t) are shown in the inset plot of Fig 1(a) 


In order to study the temporal synchronization behavior, we capture a serial 
snapshots of network state. Since the W(t) for critical transmission delay 
is the combination of a decreasing trend and periodic fluctuation (see Fig. 
1(b)), the times series ..., t m \ are chosen as W (f,) is the local maximum. 


For simplicity to compare the temporal synchronization behaviors, snapshots 
of synchronization process in the stable region, whose synchronization degree 
is indicated by W'(t), are chosen at the time series [t[, t ' 2 ,..., t' m ] satisfying 
W(ti) = W'it'j). Only snapshots of ti after 3 x 10 4 steps are considered to 
avoid the disturbance in the transient time. Total M = 356 snapshots are 


selected in each synchronization, part of which is shown in Fig. 1(b) 
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Figure 1: (Color online) (a) The logarithmic synchronizing speed of delay-coupled oscil¬ 
lates. The critical transmission delay t c = 4.2 x 10 -3 is marked by the dashed line. For 
the inner plot, it gives temporal evolutions of W(t) when r is 2.0 x 10 -3 (blue line) and 
4.3 x 10 -3 (red line), (b) Part of the snapshots of synchronization process. ■ indicates 
snapshot of synchronization on the edge of stability and • describes the corresponding one 
of synchronization on the stable region, which shows the same synchronization error. Note 
the break in x-axis. 


3. Synchronizing Path on the Edge of Stability 

According to the method in Ref. {23|, we investigate the synchronizing 
path via examining how the synchronized clusters grow in the snapshots. 
The synchronized clusters are extracted as follows. The edge linking node i 
and j is weighted by the Euclidean distance between the nodes 

»•{,., }< /,• x, - x ; (5) 


Then a threshold T is set and the edges associated to weights less than the 
threshold are preserved to form a new network for each snapshot. In this 
new networks, two connected nodes are considered as synchronized. The 
largest synchronized cluster is defined as the giant component (GC) of the 
new network. Here we choose the synchronization error of the last snapshot 
W(tM) = W\t' M ) as the threshold. The percents p(t) of preserved edges 
and sizes of GCs S(t) as a function of synchronization error W(t) for two 
synchronization processes are show in Fig. 2(a) and 2(b), respectively. Note 
that the decreasing order in x-axis of both figures is in correspondence to the 
increase of time and decrease of synchronization error. One can see that at 
the same degree of synchronization, more edges are preserved and the GCs 
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Figure 2: (Color online) (a) The percent of preserved edges, (b) The size of GC. Both 
of them imply the difference of synchronizing path between these two synchronization 
processes. Note that in both figures, the decreasing order in x-axis indicates the increase 
of time and decrease of synchronization error. 


are larger in the snapshots of synchronization with the critical transmission 
delay than those in the snapshots of synchronization on the stable region. 
This result suggests that at the same degree of synchronization the network is 
more synchronized on the edge of stability and imply a different synchronizing 
path. 

It has been observed that nodes’ degrees are tightly connected to the syn¬ 
chronization process, e.g. phase synchronization starts at high-degree nodes 


23| and explosive synchronization emergences if the natural frequencies of 

Naturally, we examine the de¬ 


oscillators depend on their degrees [2 
grees (corresponding to BA network) of nodes in the GC to give a details of 
synchronizing path. In Fig. 3(a) it shows the mean degree ( k) of nodes in 


the GC as function of its size S(t). The referred baseline is the mean degree 
k — 10 of the BA network. Concretely, when the size of GC is small, the 
mean degrees of the synchronization on the stable region are much higher 
than k = 10, which suggests that the synchronization starts at the high- 
degree nodes. As the size of GC grows, the the mean degree decreases to 
k = 10 because of the join of low-degree nodes, on the contrary, the syn¬ 
chronization on the edge of stability starts with some low-degree nodes since 
the mean degree is around k at the beginning, then the mean degree further 
decreases to smaller than k because more low-degree nodes join. It starts 
to increase till the GC is relevantly large which implies that the high-degree 
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Figure 3: (Color online) (a) The mean degree (k) of nodes in the GC, (b) The largest 
degree k max of nodes in the GC, (c) the average path length of GC, and (d) the cluster 
coefficient of GC as function of its size S(t). The dotted line indicates the mean degree, 
largest degree, average path length and cluster coefficient of network in (a), (b), (c) and 
(d), respectively. Note that in all figures, the increasing order in x-axis indicates the 
increase of time and decrease of synchronization error. 
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nodes join the GC very late. To further determine the result, we present 
the largest degree k max in the GC as a function of S(t), shown in Fig. 3(b) 
It can be seen that the largest degree K max = 121 of the BA network join 
the GC at the very beginning for transmission delay r = 2.0 x 10~ 3 , while 
it is isolated for critical transmission delay. In the BA network, high-degree 
nodes act as structural hubs and play an important role in dynamic pro¬ 
cesses such as pinning control 29] and synchronization j30(]. However, they 


lost their roles as hubs in the abnormal synchronizing path, which suggests 
the disagreement between structure and dynamics of network. As it’s unam¬ 
biguously agreed that many networks realize their functions via dynamical 
processes, the disagreement will lead to disorder of system. 

As a result of the hubs’ absence, the GCs in the synchronization on the 


edge of stability display poor small-world properties. As shown in Figs, 3(c) 


and 3(d) it is characterized by the larger average path lengths and lower 


cluster coefficients compared with their baselines of the BA network. How¬ 
ever, for the synchronization on the stable region, the average path lengths 
are shorter and cluster coefficients are higher in the normal synchronizing 
path. This variation of small-world property is also observed in the collec¬ 
tive dynamics of real complex systems, e.g. the brain, which usually suggests 
disorders in the neural activities. In an study based on fMRI of Alzheimer dis¬ 
ease (AD), clustering was significantly reduced in brain functional networks 
and the loss of small-world property is able to discriminate AD patients from 


age-matched comparison subjects with high specificity and sensitivity 31 


In another EEG study of AD, path lengths in brain functional networks con¬ 
structed by f3 -band (15 — 35 Hz) waves were significantly increased in AD 


patients 32| . Therefore, the observation of small-world properties further 


suggests the abnormal synchronizing path is potentially related to disorders 
in real complex systems. 

To investigate how the GCs diffuse to the whole network in synchroniza¬ 
tion with transmission delay, we examine whether the nodes newly joining 
the GC at each snapshot are isolated ones or in synchronized clusters at the 
previous snapshot. Assuming V[k) is the set of the nodes newly joining the 
GC at the k -th snapshot and V'(k,Ci ) is the subset of V(k) containing the 
nodes that belong to synchronized clusters labeled by c % at the (k — l)th 
snapshot. Denote the number of newly joined nodes as A S(k) = |H(/c)| and 
the maximum size of subsets as s max (k) = max* \V'(k, Cj)|. The points of 
(A S(k),s max (k)) of two synchronizing paths are shown in Fig. 4(a) and 4(b) 
respectively. Color of points indicates their number of occurrence. Small 
















Figure 4: (Color online) (a) A S(k) vs. s max of normal synchronizing path (b) A S(k) vs. 
Smax{k ) of abnormal synchronizing path (c) the growth of GC in normal synchronizing 
path, A S(k) = 20, s max {k) = 14 and k = 259. The size of GC increases from 74 to 94. (d) 
the growth of GC in abnormal synchronizing path, A S(k) = s m ax(k) = 76 and k = 301. 
The size of GC increases from 878 to 954. 
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(A S(k), s max (k )) pairs represent the nodes connect to GC as isolated ones, 
which takes most parts of (A S(k), s max (k )) pairs in normal synchronizing 
path (see Fig. 4(a)). High A S[k) and A S[k) ~ s max (k) imply the nodes join 
the GC as one synchronized cluster, which can be found in the abnormal 
synchronizing path (see Fig. 4(b)). Fig. 4(c) and 4(d) respectively show 
two examples of GC’s growth , whose A S(k) and s max (k ) are both large, in 
normal and abnormal synchronizing paths via drawing network using Pajek 
33]. In both figures, the separated clusters in the previous snapshot are indi¬ 


cated by different colors. The clustering synchronization behavior is clearly 
demonstrated in Fig. 4(d), though the communities don’t exist in the BA 


network. The result further supports the disagreement between structure 
and dynamics in abnormal synchronizing path. 


4. Conclusion 

In summary, we have shown that transmission delay can deeply influence 
the temporal synchronization behavior. Even the synchronization processes 
are all exponential convergence on the edge of stability and on the stable 
region, they display completely different synchronizing paths. In the syn¬ 
chronization with critical transmission delay, the high-degree nodes loss their 
roles as hubs in the collective dynamics. As a result, the synchronized part 
of network displays poor small-world properties in the abnormal synchroniz¬ 
ing path. In further, a detailed investigation of the growth of the largest 
synchronized clusters suggests that the clustering synchronization behaviors 
are contained in abnormal synchronizing path though the communities don’t 
exist in BA model. Our work demonstrates that the same network structure 
can display different dynamic behaviors with varying transmission delay, and 
the transition of synchronizing path leads to disagreement between structure 
and dynamics, which is probably related to the disorder of real complex 
systems, e.g. the disease of human brain. 
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